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LONG-TERM  GOALS 

The  long  term  goal  of  this  project  is  to  develop  a  better  understanding  of  oceanic  processes  in  the 
range  of  100  m  to  10  km,  in  the  so-called  submesoscale  range.  In  particular,  it  is  important  to 
explore  and  find  out  whether  and  what  type  of  submesoscale  instabilities  exist,  how  they  are 
connected  to  both  larger  scale  and  smaller  scale  motions,  and  to  what  extent  they  influence 
transport  processes  in  the  ocean.  Another  important  objective  of  this  project  is  to  test  how  well 
subgrid-scale  (SGS)  models  for  large  eddy  simulations  (LES)  work  in  the  presence  of  backward 
energy  cascade  that  may  be  characteristic  in  submesoscale  motions.  A  long  term  objective  of  this 
effort  would  be  to  improve  the  predictive  skill  of  the  Navy  numerical  models  for  submesoscale 
transport  in  the  ocean. 

OBJECTIVES 

My  main  objective  has  been  to  model  upper  ocean  mixed  layer  instabilities,  investigate  their 
behavior  and  try  to  develop  sampling  strategies  using  synthetic  drifters  and  tracers  prior  to  the 
first  LatMix  cruises.  In  addition,  modeling  help  with  data  analysis  is  planned. 

APPROACH 

The  work  is  based  on  LES  using  the  non-hydrostatic  spectral  element  model  Nek5000  as  well  as 
field  data.  Also,  software  has  been  developed  in  order  to  automatically  download  and  make 
available  MODIS  raw  SST  data  over  the  selected  region  in  the  first  LatMix  cruise. 

WORK  COMPLETED 

I  have  completed  work  on  eight  papers  (see  Publications  section)  under  the  support  of  the 
LatMix  project.  Of  these,  two  (Ozgokmen  et  ah,  2011;  Ozgokmen  and  Lischer,  2012)  contain  a 
detailed  analysis  of  mixed  layer  instability  (MLI)  with  particular  emphasis  on  sampling  strategies 
using  LES.  Two  papers  focus  on  subgrid-scale  models,  in  LES  (Berselli  et  ah,  201 1)  and  in 
ocean  models  (Marques  and  Ozgokmen,  2012).  Our  approach  in  Ozgokmen  et  al.  (2012)  is 
perhaps  the  first  truly  multi-scale  LES  aiming  to  dissect  the  interaction  between  the 
submesoscale  MLI  and  deep  mesoscale  baroclinic  instability,  and  how  their  multi-scale  signature 

1 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

30  SEP  2012  2' REPORT  TYPE 

3.  DATES  COVERED 

00-00-2012  to  00-00-2012 

4.  TITLE  AND  SUBTITLE 

A  Non-Fickian  Mixing  Model  for  Stratified  Turbulent  Flows 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROTECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

University  of  Miami, Rosenstiel  School  of  Marine  and  Atmoshperic 

Science, 4600  Rickenbacker  Causeway, Miami, FL, 33149 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBIECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

_ _ _  ABSTRACT 

18.  NUMBER  19a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  3S 

unclassified  unclassified  unclassified  Report  (SAR) 

12 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


can  be  detected  on  the  basis  of  scale-dependent  relative  dispersion  plots.  I  also  investigated  the 
observation  periods  needed  for  the  effects  of  submesoscale  and  mesoscale-induced  dispersion  to 
show  up  in  these  plots.  The  study  by  Mensa  et  al.  (2012)  is  on  the  seasonality  of  submesoscale 
features  in  the  Gulf  Stream  from  realistic,  multiple  resolution  HYCOM  simulations,  and  how  a 
combination  of  MLI  and  frontogenesis  act  to  fonn  these  features.  In  Schroeder  et  al  (2012),  we 
summarize  results  from  a  targeted  drifter  launch  at  a  coastal  frontal  zone.  In  Haza  et  al.  (2012), 
we  evaluate  two  different  classes  of  subgrid-scale  parameterizations  for  Lagrangian  dispersion  in 
HYCOM  for  the  North  Atlantic  in  order  to  compensate  for  the  effect  of  enhanced  dispersion  over 
submesoscale  regime  that  is  typically  missing  in  mesoscale  eddy  resolving  models. 

RESULTS 

A  subset  of  results  are  presented  below,  one  from  a  computational  LES  study  and  another  from  a 
field  experiment  on  lateral  mixing.  I  have  presented  a  number  of  new  computations  during  the 
LatMix  meeting  in  June  2012,  which  have  not  been  prepared  for  publication,  and  therefore,  not 
included  below. 

1)  Deduction  of  submesoscale  and  mesoscale  interactions  from  dispersion: 

This  study  is  motivated  by  a  recent  observation  of  the  emergence  of  submesoscale  features  along 
the  rim  of  a  cyclonic  Gulf  Stream  ring  that  was  noted  during  the  2011  Summer  LatMix  cruise 
(Fig.  1).  These  submesoscale  features  developed  within  days  after  a  cyclonic  (cold  core)  ring 
entrained  warm  water  from  the  Gulf  Stream,  but  apparently  completely  dissipated  within  three 
weeks  of  their  emergence.  Submesoscale  features  have  been  notoriously  difficult  to  observe  in 
situ,  while  remote  sensing  instruments  provide  occasional  glimpses  into  their  structure  and  rapid 
evolution.  The  implication  of  the  multi-scale  interaction  in  this  particular  case  is  that 
submesoscale  flows  appear  sporadically  because  special  circumstances  are  needed  for  their 
fonnation,  and  dissipate  quickly  under  the  action  of  mesoscale  stirring. 

During  the  time  period  when  both  the  mesoscale  eddy  and  submesoscale  features  coexist,  the 
question  arises  as  to  which  scales  control  transport.  The  primary  objective  of  this  study  is  to 
investigate  the  validity  of  two  contrasting  hypotheses:  Hypothesis-I  where  mesoscale  eddies 
exert  non-local  control  over  small  scale  transport  and  Hypothesis-II  where  dispersion  and 
transport  at  the  submesoscales  is  governed  by  eddying  motions  at  local  scales.  The  implication  of 
Hypothesis-I  is  that  mesoscale  eddy-resolving  data-assimilating  (current  generation)  OGCMs 
would  provide  adequate  predictions  in  applied  important  problems  such  as  transport  of  coastal 
pollutants  and  oil  spills,  while  Hypothesis-II  would  indicate  that  OGCMs  would  significantly 
underestimate  the  dispersion  of  a  passive  pollutant  until  the  patch  size  grows  to  the  scale  of  the 
local  radius  of  defonnation.  Another  complementary  objective  of  this  study  is  to  provide  insight 
into  the  question  of  whether  and  how  instabilities  below  the  mixed  layer  influence  surface 
transport. 

LES  is  used  to  assess  transport  properties  of  multi-scale  ocean  flows.  In  particular,  the  initial 
density  perturbation  field  is  specified  in  a  way  that  MLI  eddies  are  approximately  an  order  of 
magnitude  smaller  than  deeper  eddies.  A  parallel  experiment  is  configured  for  comparison  in 
which  the  stratification  below  the  mixed  layer  has  no  available  potential  energy,  and  thereby  no 
independent  turbulent  coherent  structures  from  the  surface  motion  can  form. 
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In  particular,  we  compare  scale-dependent  measures  of  Lagrangian  relative  dispersion  and  the 
evolution  of  passive  tracer  releases  in  numerical  experiments  containing  only  submesoscale  MLI, 
and  MLI  modified  by  deeper,  baroclinic  mesoscale  disturbances. 

Visualization  of  turbulent  coherent  structures  through  3D  finite-time  Lyapunov  exponents 
(FTLEs)  and  analysis  of  surface  kinetic  energy  wave-number  spectra  show  that  smaller  scale 
instabilities  in  the  mixed  layer  rapidly  lose  their  coherence  in  the  presence  of  larger  scale 
straining  induced  by  mesoscale  motions  (Fig.  2).  Tracer-based  eddy  diffusivity  computed  to 
quantify  the  turbulent  exchange  across  the  mixed  layer  front  yields  approximately  0.25  m“/s 
uniformly  within  the  mixed  layer  depth  in  the  case  of  MLI  only,  but  the  frontal  turbulent 
exchange  diffusivity  increases  to  5  m  /s  as  the  flow  transitions  to  mesoscales  (Fig.  3,  upper 
panel).  This  result  is  likely  to  be  related  to  shallow  mixed  layer  depths  representative  of  summer 
conditions  in  the  North  Atlantic.  Winter  mixed  layers  are  typically  an  order  of  magnitude 
deeper,  possibly  creating  large  enough  mixed  layer  eddies  that  can  survive  the  larger  scale 
straining  for  long  periods. 

The  scale-dependent  FSLE  /.(8)  curves  for  both  cases  are  computed  by  releasing  7803  synthetic 
Lagrangian  particles  within  a  5  km  by  5  km  square  over  the  front  at  three  levels  within  the 
mixed  layer.  The  numerical  experiment  in  which  only  the  MLI  is  present  leads  to  results 
consistent  with  Hypothesis-II,  namely  a  single  FSLE  plateau  emerges  at  separation  scales 
smaller  than  the  MLI  eddies.  Given  the  small  length-scales  and  energy  of  the  MLI  eddies, 
however,  the  inverse  cascade  produces  a  relatively  weak  mesoscale  flow  over  the  limited  time 
period  examined.  The  separation  times  increase  approximately  linearly  with  separation  scale  for 
all  scales  greater  than  approximately  300  m. 

The  primary  new  advance  made  by  this  study  is  that  during  the  period  when  both  MLI  and 
deeper  instabilities  are  present,  X  (8)  at  0  and  10  m  depths  show  two  distinct  plateau  clearly 
associated  with  the  disparate  instability  scales  (Fig.  3,  lower  panel).  There  is  a  complete  collapse 
of  these  curves  with  the  one  calculated  from  particle  releases  at  25  m  for  separation  distances 
larger  than  the  scale  of  MLI  eddies  (8  >  1000  m),  indicating  that  larger  scale  motions  are 
controlled  by  deep  instabilities.  While  we  do  not  directly  compare  flows  with  and  without 
submesoscale  features,  the  strong  effect  of  such  small-scale  features  on  surface  transport  can  be 
inferred  from  the  significant  differences  in  the  dispersion  characteristics  observed  above  and 
below  the  mixed  layer  depth  in  Exp  II. 

These  results  indicate  the  distinct  possibility  of  Hypothesis-II  holding,  at  least  episodically,  in 
the  upper  ocean.  For  computational  reasons,  we  have  examined  a  situation  with  two  disparate 
length  scales  imposed.  One  could,  however,  easily  imagine  multiple  mixed  layer  fronts  driven 
by  a  range  of  lateral  buoyancy  gradients  forming  eddies  of  various  scales  with  their  mutual  inter¬ 
action  creating  a  nearly  continuous  spectrum  of  submesoscale  features.  We  have  put  forward 
parameterizations  to  compensate  for  the  dispersion  defect  between  Hypothesis-I  and  Hypothesis- 
II  for  the  submesoscales  in  OGCMs  have  been  put  forward  recently  (Haza  et  ah,  2012). 

Given  that  simultaneous  release  thousands  of  drifters  is  currently  not  feasible  in  oceanic 
experiments,  but  numbers  on  the  order  of  one  hundred  can  be  readily  considered  in  dedicated 
studies,  the  sensitivity  of  X  (8)  is  explored  to  the  number  of  drifters  and  duration  of  the 
trajectories  so  that  both  plateau  can  be  captured.  The  results  indicate  that  the  submesoscale 
plateau  emerges  readily  with  even  32  drifters,  while  timing  and  duration  are  important 
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considerations  to  reveal  both  plateau.  Timing  is  important  in  that  submesoscales  do  not  coexist 
with  mesoscales  perpetually  while  the  drifters  must  sample  the  flow  field  long  enough  to 
separate  under  slower  evolving  mesoscale  motions.  The  results  imply  that  a  field  experiment 
would  require  fast  sampling  drifters  (on  the  order  of  minutes)  and  at  least  one-month  long 
trajectories  to  span  the  scale  separation. 

2)  Targeted  sampling  of  submesoscale  dispersion  at  a  coastal  frontal  zone: 

Relative  dispersion  at  oceanic  submesoscales  is  still  largely  unknown,  and  has  been 
experimentally  investigated  so  far  mostly  in  the  open  ocean  and  for  scales  greater  than  1  km 
(LaCasce  and  Ohhnann,  2003;  Koszalka  et  ah,  2009,  Lumpkin  et  ah,  2010,  Schroeder  et  ah, 

2011,  Berti  et  ah,  2011).  With  collaborators  in  Italy  and  France,  we  have  conducted  two 
experiments  in  the  summer  of  2010,  which  directly  target  relative  dispersion  induced  by  a 
persistent  coastal  front  in  the  submesoscale  range  of  100  m  to  1000  m.  A  novel  sampling 
strategy  has  been  used  to  provide  insight  into  the  quantification  of  initial  tracer  dispersion  under 
submesoscale  motions. 

The  area  of  interest  is  the  Eastern  Ligurian  Sea  (in  the  North-Western  Mediterranean  Sea,  Fig.  4 
upper  right).  The  region  is  characterized  by  depths  shallower  than  200  m  and  is  influenced  by 
the  Amo  river  discharge  along  with  some  minor  rivers,  creating  favorable  conditions  for  the 
fonnation  of  buoyancy  fronts  within  the  coastal  area.  The  Eastern  Ligurian  Sea  is  dynamically 
different  from  the  central  and  deeper  part  of  the  Ligurian  basin,  which  is  characterized  by  a 
strong  and  coherent  cyclonic  circulation  (the  Liguro-Provencal  current),  with  a  typical 
defonnation  radius  of  10  km  or  more.  Relative  dispersion  in  the  Liguro-Provencal  current  was 
previously  studied  within  the  framework  of  two  field  experiments  (Marine  Rapid  Environmental 
Assessment,  MREA07/08),  where  drifters  were  released  in  clusters  with  a  separation  distance  of 
1  km  (Schroeder  et  ah,  2011).  MREA07/08  results  showed  the  presence  of  a  non-local  dispersion 
regime  for  the  submesocale  range  of  1  to  10  km. 

The  two  experiments  reported  here  took  place  in  the  Eastern  Ligurian  Sea  during  the  summer  of 
2010  and  were  conducted  by  an  extensive  international  partnership.  The  Ligurian  Dispersion 
Experiment  (hereafter  LIDEX10)  was  led  by  the  Italian  Consiglio  Nazionale  delle  Ricerche 
during  July  2010  with  the  RV  Maria  Grazia,  while  the  Recognized  Environmental  Picture  2010 
experiment  (hereafter  REP  10)  was  led  by  the  NATO  Undersea  Research  Center  with  the  NRV 
Alliance  during  August-September  2010.  During  both  experiments,  drifter  deployment 
locations  were  selected  with  the  help  of  satellite  images,  especially  the  chlorophyll  concentration 
(Fig.  4,  lower  left)  to  identify  the  frontal  area  separating  rich/turbid  coastal  waters  from  the 
offshore  waters.  In  the  case  of  LIDEX10,  in  situ  salinity  data  were  collected  by  performing  CTD 
casts  every  0.5  to  1  km  across  the  front  (from  the  surface  to  a  depth  of  50  m),  while  a  Scanfish,  a 
towed  undulating  vehicle,  was  used  to  precisely  delineate  the  front  during  REP  10 
(Fig.  4,  lower  right). 

The  data  set  used  to  compute  relative  dispersion,  the  primary  metric  of  interest,  derive  from 
CODE  drifters  deployed  with  a  multi-scale  approach.  The  deployment  kernel  consisted  of 
clusters  of  nine  drifters,  composed  of  three  triplets,  initially  separated  by  distances  ranging 
between  50  m  and  600  m.  This  strategy  was  chosen  to  assure  that  a  sufficiently  high  number  of 
pairs  at  different  separation  distances  would  be  available  for  the  computation  of  dispersion 
metrics.  Nineteen  drifters  were  launched  during  LIDEX10  and  eleven  drifters  during  REP  10. 
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To  support  the  interpretation  of  the  ongoing  dynamics,  the  ocean  model  ROMS  was  employed  to 
simulate  the  Ligurian  Sea  circulation  during  the  LID  EX  10  and  REP  10  trials.  The  horizontal 
resolution  of  this  model  is  1 .7  km,  with  32  vertical  sigma-coordinate  levels  (see  Fig.  4  for  the 
Ligurian  Sea  ROMS  domain).  The  model  was  initialized  in  May  2010  with  the  Mediterranean 
Forecasting  System  (MFS)  model  analysis.  Open  boundary  conditions,  provided  by  MFS,  were 
applied  to  tracers  and  baroclinic  velocity  with  radiation  and  nudging,  while  free  surface  and 
depth-integrated  velocity  boundary  conditions  were  applied  following.  The  high-resolution,  non¬ 
hydrostatic,  atmospheric  model  COSMO-ME  of  the  Italian  Air  Force  National  Meteorological 
Center  (Centro  Nazionale  per  la  Meteorologia  e  Climatologia  Aeronautica)  provided  the  surface 
forcing.  Three  major  rivers  are  included  as  lateral  boundary  conditions:  Arno,  Serchio  and  Magra 
(monthly  climatologies).  No  data  assimilation  was  employed. 

To  study  the  dispersion  regimes  from  the  ocean  model,  synthetic  drifters  were  released  in  the 
model,  following  two  strategies:  (a)  Synthetic  drifters  are  launched  corresponding  to  the  actual 
time  and  location  of  real  deployments  LIDEX10  and  REP  10.  (b)  In  order  to  ensure  statistical 
reliability  of  the  results,  large  ensembles  of  synthetic  trajectories  are  also  used  for  both 
LID  EX  10  and  REP  10  experiments.  The  ensembles  consisted  of  135  particles  launched  1 1  times 
at  12  hours  intervals  around  the  time  of  actual  in-situ  drifter  deployments.  The  spatial 
distribution  of  each  launch  consisted  of  29  grid  points  centered  around  the  in-situ  drifter 
launching  points,  with  a  cross-configuration  corresponding  to  1  center  particle  at  each  grid 
point  and  4  satellite  ones,  displaced  by  500  m  along  the  four  cardinal  directions.  Drifter 
trajectories  were  integrated  on-line  in  2D  at  0.5  m  depth  in  order  to  mimic  the  response  of  the 
real  drifters,  and  trajectories  sampled  at  30  minute  intervals  are  used  to  compute  relative 
dispersion  statistics. 

The  scale  dependent  FSLE  highlights  relative  dispersion  over  the  submesoscales.  Results  from 
LID  EX  10  and  REP  10  are  presented  in  Fig.  5  and  clearly  demonstrate  a  dispersion  regime  that  is 
scale-dependent  (local)  at  all  separation  scales.  Results  from  both  LIDEX10  and  REP  10  data  are 
consistent  with  one  another.  MREA07/08  results  ranging  down  to  8  ~  100  m  are  also  shown,  on 
the  basis  of  chance  pairs.  The  contrast  with  respect  to  MREA07/08  is  quite  significant  over  the 
submesoscales  in  that  there  is  a  persistent  increase  in  the  FSLE  with  a  Richardson's  power  law  of 
X  ~  8'  over  the  submesoscales  (8  <  10  km)  ending  up  with  a  limiting  value  of  Amax~  20  days' 
at  the  smallest  separations  of  8<100  m  attained  during  LIDEX10  and  REP10.  This  value  of  the 
limiting  FSLE  is  about  an  order  of  magnitude  greater  than  that  obtained  in  the  Ligurian  Sea 
circulation  during  MREA07/08.  ROMS  simulations  clearly  lack  any  dispersive  turbulent 
motions  below  the  mesoscale  separations  of  8  <  10  km.  There  is  however  very  good  agreement 
between  results  from  launches  with  strategy-(a)  and  strategy-(b)  for  LID  EX  10  over  the 
submesoscales  (same  for  REP  10,  not  shown).  Note  that  all  curves  converge  over  the  mesoscale 
regime  of  8>10  km  and  show  ballistic  scaling  X  ~  8': ,  associated  with  horizontal  shear  zones 
created  by  gyres  and  boundary  currents  in  the  Ligurian  Sea. 

It  is  concluded  that  tightly-spaced  clusters  of  drifters  constitute  an  efficient  and  inexpensive 
observational  technique  to  obtain  relative  dispersion  statistics  from  rapidly-evolving 
submesoscale  flows  with  minimal  aliasing.  It  is  found  on  the  basis  of  both  relative  dispersion  and 
and  scale-dependent  FSLE  that  time  scales  of  t<15  hours  and  space  scales  of  8<1  km,  that  are 
characteristic  of  submesoscale  phenomena,  exhibit  enhanced  relative  dispersion  compared  to 
larger  scales  of  motion.  In  particular,  values  were  obtained  reaching  20  days'1  for  8  <  100 
m,  which  is  an  order  of  magnitude  greater  than  those  obtained  for  the  off-shore  mesoscale  flows. 
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Computations  with  realistically-configured  1 .7  km  horizontal  resolution  ROMS  were  unable  to 
capture  this  enhanced  dispersion  regime,  indicating  the  need  for  subgrid-scale  models  for 
submesoscale  transport. 


IMPACT/APPLICATIONS 

The  scales  considered  in  this  project  represent  the  range  of  scale  of  navy  operations  and  thus 

anomalous  currents  and  perturbations  in  the  acoustic  and  optical  environment  that  can  affect  a 

variety  of  navy  operations.  Understanding  the  motion  in  this  range  of  scales  is  therefore  critical 

to  help  improve  the  predictive  capability  of  the  existing  Navy  models. 
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Fig.  1:  An  example  of  submesoscale  instabilities  from  MODIS  SST  images  in  the  Gulf 
Stream  region  during  the  2011  Summer  LATMIX cruise.  (Upper  left)  Entrainment  of 
a  warm  filament  around  a  cyclonic  eddy,  t  =  0  (May  1,  2011).  (Upper  right)  The  onset 
what  appears  to  be  mixed  layer  instabilities  att  =  1  day.  (Lower  left)  Formation  of 
sinusoidal  perturbations  along  the  rim  of  the  cyclonic  eddy  att  =  2  days.  (Lower 
right)  The  state  of  the  eddy  at  t  =  22  days. 
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Fig.  2:  Density  perturbation  fields  from  LES  and  3D  FTLEs  (yellow)  in  the  case  of 
MLI  only  (upper  panel)  and  when  also  deep  baroclinic  instability  is  also  present  (lower 
panel).  Note  the  disappearance  of  MLI  after  deep  baroclinic  eddies  set  in.  The 

animations  are  available  from: 

http://www.rsmas.miami.edu/personal/tamay/3D/mlisl3.mov  and 
http://www.rsmas.miamledu/personal/tamay/3D/mlisll.mov 
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Fig.  3:  (Upper  panel)  Time  evolution  of  the  second  moment  of passive  scalar  at  the 
surface  and  at  the  base  of  the  mixed  layer  in  simulations  containing  only  the  MLI 
(Exp-I,  solid  lines)  and  also  the  deep  baroclinic  eddies  (Exp-II,  dashed  lines).  (Lower 
panel)  Scale-dependent  FSLE  from  particle  releases  at  the  surface,  middle  and  base  of 
the  mixed  layer.  Note  in  Exp-II  the  the  double  plateau  within  the  mixed  layer,  one 
associated  with  mesoscale  eddies  and  the  other  with  submesoscale  MLI  eddies. 


10 


□ 

□ 


Domain  of  ROMS 

Ocean  colour  domain  (below) 
REP10  cluster  deployments 


mg/mA3 


43.7 


10  1 


43  65 


10.15 


Latitude  N 


43.6  10.05  Longitude  E 


■  -37.9 


37.8 

37.7 

37.6 

37.5 


Fig.  4:  (Upper  right)  Region  of  interest,  numerical  model  domain  and  and  satelitte 
image  frame.  (Lower  left)  Ocean  color  image  (chll  retrieved  from  MERIS  09:33  UTC)  of 
21st  August  2010  (before  the  deployment  of  the  REP10  drifters)  with  the  REP10  cluster 
deployment  positions  and  the  track  of  the  ScanFish  MKII  used  for  in  situ  data  collection. 
(Lower  right)  Vertical  salinity  section  along  the  ScanFish  MKLI  track  with  the  REP10 
cluster  deployment  positions  across  the  frontal  region. 
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Fig.  5:  Scale  dependent  FSLE  for  LIDEX10  (solid  blue),  REP10  (solid  green), 
and  MREA07/08  (solid  red).  Modeling  results  from  synthetic  drifters  in  ROMS  simulations 
are  shown  as  follows:  19  drifters  launched  as  in  the  LIDEX10  (blue  with  “+”, 
as  in  real  in-situ  experiment),  1485  drifters  during  LIDEX10  (dashed  blue,  high  number  of 
synthetic  drifters  in  LIDEX10)  and  1485  drifters  during  REP10  (dashed  green).  Ballistic 

(A~S  )  and  Richardson  (A~5  )  regimes  are  indicated. 
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